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Abstract. We consider a linear size-structured population model with diffusion in 
the size-space. Individuals are recruited into the population at arbitrary sizes. The 
model is equipped with generalized Wentzcll-Robin (or dynamic) boundary condi- 
tions. This allows modelling of "adhesion" at extremely small or large sizes. We 
establish existence and positivity of solutions by showing that solutions are governed 
by a positive quasicontractivc semigroup of linear operators on the biologically rel- 
evant state space. This is carried out via establishing dissipativity of a suitably 
perturbed semigroup generator. We also show that solutions of the model exhibit 
balanced exponential growth, that is our model admits a finite dimensional global 
attractor. In case of strictly positive fertility we are able to establish that solutions 
in fact exhibit asynchronous exponential growth. 



1. Introduction 

Starting with the seminal work of Gurtin and MacCamy |17j . a significant amount 
of interest has been devoted to the analysis of mathematical models arising in struc- 
tured population dynamics (see e.g. [31] US] for further references). Such models often 
assume spatial homogeneity of the population in a given habitat and only focus on the 
dynamics of the population arising from differences between individuals with respect to 
some physiological structure. In this context, reproduction, death and growth character- 
ize individual behavior which may be affected by competition, for example for available 
resources. 

In a recent paper, Hadeler |19j introduced size-structured population models with dif- 
fusion in the size-space. The biological motivation is that diffusion allows for "stochastic 
noise" to be incorporated in the model equations in a deterministic fashion. The main 
question addressed in |19j is what type of boundary conditions are necessary for a bio- 
logically plausible and mathematically sound model. In this context some special cases 
of a general Robin boundary condition were considered. 
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In this paper we introduce the following linear size-structured population model 

"i( s > t) + (7(s)u(s, t)) s = (d{s)u s (s, t)) s - fj,(s)u(s, t)+ P(s, y)u{y, t) Ay, s e (0, m), 

Jo 

(1.1) 

[(d(s)u,(«, *)) J s=0 - & ou s (0, t) + c u(0, t) = 0, (1.2) 
[(d(s)u s (s, t)) s ] s=m + b m u s (m, t) + c m u(m, t) = 0, (1.3) 



with a suitable initial condition. The function u = u(s, t) denotes the density of indi- 
viduals of size, or other developmental stage, s at time t. The non-local integral term 
in ([l.ip represents the recruitment of individuals into the population. Individuals may 
have different sizes at birth and f3(s,y) denotes the rate at which individuals of size y 
"produce" individuals of size s. fx denotes the size-specific mortality rate while 7 denotes 
the growth rate, d stands for the size-specific diffusion coefficient, which we assume to 
be strictly positive. 60 and b rn are positive numbers, while cq and c m are non-negative. 
We will discuss special values for these constants in Equation (|1.7[) below. 

Equation (|1.1[) describes the evolution of a "proper" size-structured population, in 
contrast to the one where it is assumed that all newborns enter the population at a 
minimal size. In that case the linear model can be rewritten as an age-structured model 
using the Grabosch-Hcijmans transform, see e.g. [21]. We also refer the interested reader 
to [3j [3 [10] where size-structured models with distributed recruitment processes were 
investigated. The boundary conditions (|1.2p - (|1.3[) are the so called generalized Wentzell- 
Robin or dynamic boundary conditions. These "unusual" boundary conditions were 
investigated recently for models describing physical processes such as diffusion and wave 
propagation, see e.g. [TTJ [T^] [TB] . Briefly, they are used to model processes where particles 
reaching the boundary of a domain can be either reflected from the boundary or they can 
be absorbed and then released after some time. Our goal here is to introduce this rather 
general type of boundary condition in the context of models describing the evolution of 
biological populations, with particular focus on positivity and asymptotic behavior of 
solutions. 

The first works introducing boundary conditions that involve second order derivatives 
for parabolic or elliptic differential operators go back to the 1950s, see the papers by 
Feller [T3J [T3] and Wentzell (also transliterated as "VentceF", Aleksandr D. Ventceh,) 
[2"31 l2~i] . These first studies were purely motivated from the mathematical point of view. 
The original question was to identify the set of all possible boundary conditions that 
give rise for a parabolic differential operator to generate a contraction semigroup on an 
appropriate state space. The abstract mathematical analysis gave a clue where physical 
intuition failed previously. 

The boundary conditions (|1.2[) - (|1.3[) arc in general form, and we shall now specify the 
constants bo, b m and Co, c m to give a biological explanation for the boundary conditions. 
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Integration of Equation from to m yields for U(t) = f™ u(s, t) ds 

4- U(t) = 7(0)u(0, t) - 7 (mWm, t) + d(m)u s (m, t) - d(0W0, t) 
at 

+ / / P(s,y)u(y,t)dyds - / fi(s)u(s,t)ds 
Jo Jo Jo 

=: B(t) - D(t) + 7(0)u(0, i) - j(m)u{m, t) + d(m)u s (m, t) - d(0)u s (0, i), 

(1.4) 

where i? and -D denote the combined birth and death processes, respectively. We note 
that formally, by replacing the diffusion term by its counterpart from Equation (jl.lj) . the 
boundary conditions (|1.2|) - (|1.3I) can be cast in the dynamic form 



u t (0, t) = u(0, t)(-y (0) - /x(0) - c ) + u s (0, t)(b - 7 (0)) + / /3(0, y)u(y, t) dy 



u t (m, t) = u(m, *)(— 7'(m) - j«(m) - c m ) + w s (m, i) (-6 m - j(m)) + / /3(m, y)u(y, t) dy 



o 

(1.5) 

These are the governing equations for individuals of minimum and maximum sizes, re- 
spectively. It is natural to assume that in the absence of mortality and recruitment, 
i.e. when B{-) = D(-) = 0, the total population size remains constant at every time t. 
Mathematically, this amounts to the condition 

0= j t U{t) + ut(0,t)+u t {m,t) 

= u(0, t) ( 7 (0) - 7 '(0) - c ) + u(m, t) (- 7 (m) - c m - 7 '(m)) ( L6 ) 
+ u,(0, t) (b - d(0) - 7 (0)) + u a (m, t) (d(m) - j(m) - 6 m ) . 

In order to guarantee conservation of total population in the absence of birth and death 
processes, we make the following assumptions 

c = 7 (0) - 7 '(0), c m = - 7 (m) - 7 '(m), 6 = d(0) + 7 (0), b m = d(m) - j(m). (1.7) 

We note that condition (jl.7p together with the assumption that b is positive and c is non- 
negative impose a restriction on the growth rate 7. The boundary conditions (|1.2|) - (|1.3|) 
can now be written as 

d(0)u ss (0, t) + (d'(0) - d(0) - 7 (0)) u a (0, *) + (7(0) - 7 '(0)) u(0, t) = 0, 

d(m)u ss (m, t) + (d (m) + <f(m) — 7(m)) u s (m, t) + (—7(771) — j'(m)) u(m, t) = 0. 

The boundary conditions have now the following natural interpretation. Individuals in 
the population whose size becomes very small or very large can be absorbed in the states 
s = or s = m and become inactive, while after some time they become active again 
due to diffusion. We further impose the following assumptions on the model ingredients 

/i G C([0,to]), /3 G C([0,m] x [0,m]), 0,(j,>Q, 7,d G C\[0,m\), d>0. 
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In this note first we establish existence and positivity of solutions of model (|l.ll) - (|1.3p . 
This existence proof follows similar arguments developed in [TTJ [T2]. The significant 
difference is that solutions to our model are not necessarily governed by a contraction 
semigroup, hence as in [El |5] we need to rescale the semigroup to obtain the dissipativity 
estimate. As a result of the dissipativity calculation, we show that the resolvent operator 
of the semigroup generator is positive. This was not established in [12] . 

In Section[3]we investigate the asymptotic behavior of solutions. First we establish that 
solutions of the model equations exhibit balanced exponential growth, in general. This is 
an interesting phenomenon, often observed for linear structured population models, see 
e.g. [HI HOI ■ I n some sense, it is a stability result and characterizes the global asymptotic 
behavior of solutions to the model. Then, assuming that fertility is strictly positive, we 
are able to show that after a rescaling by and exponential factor solutions actually tend 
to a fixed size-distribution. This is shown via establishing irreducibility of the governing 
semigroup. 



In this section we are going to establish the existence of a positive quasicontractive 
semigroup of operators which governs the evolution of solutions of (|1.1[) - (|1.3[) . For basic 
definitions and results used throughout this paper we refer the reader to [TJ [H [6] . Let 

X = (L 1 (0,m) ©R 2 , || • \\x), where for (x,xo,x m ) £ X the norm is given by 



for some positive constants c\ and C2 that we will specify later. Then X is a Banach lat- 
tice. Wc identify a function u £ C[0,m] with its restriction triple (u|(o,m)> u(0) 7 u(m)) £ 
X. With this identification, the set C 2 [0, m] is dense in X with respect to the 1 1 • |^-norm. 
Let 



D(A) = {u e C 2 [0, m] i^iiel^m), (d(s)u' (s))'| s=0 - b o u'(0) + cou(0) = 0, 
(d(s)u'(s))'\ 8=m + b m u(m) + c m u(m) = 0} , 



2. Existence and Positivity of Solutions 



||(a;,3;o,x m )||;t = ||a;|| L i +ci|x | + c 2 \x m \, 



(2.8) 



where 



*u(s) = (d(s)u'(s))' - (j(s)u(s)Y - /J,(s)u(s) + 




The operator A with domain D (A) is then defined by 



Au 



(b - 7 (0))£u 

{-b m -"t{m))£ s u 




where we have set 



po = /i(0) + c + 7'(0), p m = K m ) + Cm + l'{rn), 
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for short. We use [6j Theorem 3.15] (see also [2j Section A-II, Theorem 2.11]) which 
characterizes generators of contractive semigroups via dissipativity. We establish the 
existence of a quasicontractive semigroup for the general boundary condition (jl.2p ~ 
(|1.3[) . Recall that a strongly continuous semigroup {T(t)} t>0 is called quasicontractive 
if |T(i)| < e w * for some u G R and it is called contractive if the choice u> < is possible. 
A linear operator A with domain D{A) is dissipative, if for all x G D(A) and A > one 
has 1 1 (AX- A)x\\ > X\\x\\. 

Theorem 2.1. Assume that 

d(0) dim) , x 

ci = 7 V ' > 0, and c 2 = t-4 , , >0. 2.9 
oo - 7(0) 7( m ) + »m 

T/ien £/ie closure of the operator A is the infinitesimal generator of a positive quasicon- 
tractive semigroup of bounded linear operators on the state space X (where the weights 
in equation (|2.8[) are chosen accordingly). 

Proof. We introduce the modified operator A with f3 — in the definition of A. For 
A > and h G X, u 6 £*(^4) we consider the equation 

u-\(l-ul^u = h, on [0,m]. (2.10) 

That is 

h(s) = u(s) — A (^(d(s)u' (s))' — ('j(s)u(s))' — fj,(s)u(s) — Lou(s)j , s G (0, m), 

(2.11) 

A- 1 /i(0) = u'(0)(7(0) - 6o) + u(0) (A" 1 + 7 '(0) + /i(0) + c + «) , (2.12) 
A -1 /i(m) = it'(m)( 7 (m) + b m ) + u(m) (A -1 + i(m) + /u(m) + c m + w) . (2.13) 



Next wc multiply Equation (|2.11l) by Xu+( s )j where \u+ denotes the characteristic func- 
tion of u + , and integrate from to m. The boundary d[u > 0] consists of two parts, 
Ti = d[u > 0] n (0, m) and r 2 = [u > 0] n {0,™}. The term d(s)u'(s) gives negative 
contributions on Ti since the outer normal derivative of u is negative. The term 7(s)u(s) 
gives no contributions on Ti since u = there. Hence we obtain 

||" + ||i <X(sgn{u + (m))d(m)u'{m) - sgn(u+(0))d(0)u'(0)) 
— A (sgn(u + (m))7(m)w(m) — sgn(M + (0))7(0)u(0)) 

-X Xu+(s)Ks)u(s)ds - Xu\\u + \\i + / Xu+(s)h(s)ds. 
Jo Jo 
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This, combined with equations f|2.12[) - (|2.13[) . yields 

||u + ||i + sgn(u+ (m))u(m) ( (1 + A ( 7 '(m) + /x(m) + c rn + u)) + A 7 (m) j 

+ sgn( U +(0)) W (0) ( ^3^ (1 + A (7'(0) + M(0) + co + w)) - A 7 (0) 

< -A / x„+(s)^(s)w(s) ds - Aw||u + ||i + / Xu+(s)Ms)ds 
Jo Jo 

dijri) , _,„. + , n , u ,, n , d(0) 



+ sgn^^Mm);^-^ + sgn^O))/^ _ 7(Q) • 

(2.14) 

Similarly, multiplying Equation (|2.11| by — x u - ( s ) an d integrating from to m, we 
obtain 

- sgn(u~ (m))u(m) ( (1 + A ( 7 '(to) + /x(m) + c m + u)) + A 7 (m) 

\~f{m) + b m 

B gn(«-(0))tt(0) ^ o ^ (0) (1 + A ( 7 '(0) + MO) + co + «,)) - A 7 (0) 
<A/ Xu-(s)^(s)w(s)ds - Aw||m _ ||i - / x«-(s)ft(s)ds 



(2.15) 



sgn(u-(m))fe(m) ,^ W \ - sga(u~ (0))h(0)- ^ 



7 (m) + 6 m b - 7 (0) ' 

Adding ||PH) and (|2~T5l) yields 

ll«l|i + l«(o)| (^3^0) (1 + A (7 ' (0) + M0) + c ° + w)) ~ A7(0) 

/ dim) 

+ \u(m)\ - \ ' (1 + A ( 7 '(m) + /i(m) + c m + w)) + A 7 (m) 
V 7 (m) + 6 m 

< jT (x„+ (*) " Xu- (a)) h(s) ds + h(0) bo d }^ {Q) (sgn(u+(0)) - sgn(«-(0))) - MW\i 

dim) _ f m 

+ h(m) f—r—r (sgn(u + (m)) — sgn(u~(m))) - A / (x„+(s) - X„- (s)) fJ,{s)u(s) ds. 
l(m) + o m Jo 

(2.16) 

Assuming that condition (|2.9I) holds true the left hand side of inequality (|2 . 16[) can be 
estimated below, for u € M large enough, by 

IHIi + l«(o)lr^ + l«(m)| : d{m) 



6o-7(0) /l 7(m) + 6 t , 
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Similarly, the right hand side of inequality (|2.16[) can be estimated above by 

ll% + l*(o)lr^ + K»OI- d(m) 



bo - 7(0) 7(771) + b m ' 

Hence, if condition (|2.9[) is satisfied we have the dissipativity estimate 

IMI* < \\h\\ x = ||u-a(a-wx)«|| 

for the operator A — loI. 

For the range condition we need to show that whenever h is in a dense subset of X 
then the solution u of Equation (|2.10|) belongs to the domain of A. Since we assumed that 
d > 0, i.e. we have true Wentzell boundary conditions at both endpoints of the domain, 
the required regularity of the solution u € C 2 [0, m] and hence the range condition follows 
from [15j Theorem 6.31]. Thus the closure of A — uil is a generator of a contractive 
semigroup by [5J Theorem 3.15]. 

Next we observe that for h <E X + every term on the right hand side of Equation ()2,15[) 
is non-positive, while every term on the left hand side of Equation (|2.15[) is non-negative. 
The inequality can only hold for h 6 X + if = u(0) = u(m) = 0. This proves that 

the resolvent operator R(X, A — uil) = (XI — (A — ul))^ 1 is positive for A large enough, 
hence the closure of A — loI generates a positive semigroup. A simple perturbation result 
yields that the closure of A is a generator of a positive quasicontractive semigroup. 

Finally we note that, the operator A — A is positive, bounded and linear, hence it 
generates a positive semigroup {<S(f)} t>0 which satisfies 

\\S(t)\\<e tB , where B = WPW^. 

The proof of the Theorem is now completed on the grounds of the Trotter product 
formula, see e.g. O Corollary 5.8 Ch. 3]. □ 

Remark 2.2 We note that the biologically natural assumptions in (|1.7[) that guarantee 
conservation of the total population in the absence of mortality and recruitment imply 
that c\ — C2 = 1. Hence the mathematical calculations coincide with the biological 
intuition. It can be shown that in the absence of mortality and recruitment the operator 
A generates a contraction semigroup on the state space (X, || ■ \\x), in fact ||T(£)||;e = 1 
for all t > 0. 



3. Asymptotic behavior 

In this section we investigate in the framework of semigroup theory the asymptotic 
behavior of solutions of model (|1.1[) - (|1.3[) . Since our model is a linear one, we may expect 
that solutions either grow or decay exponentially, unless they are at a steady state. For 
simple structured population models in fact it is often observed (see e.g. [51 [51 [TU1 fT5] ). 
that solutions grow exponentially and tend to a finite-dimensional (rcsp. one-dimensional) 
global attractor. This phenomenon is called balanced (resp. asynchronous) exponential 
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growth. The rate of exponential growth is called the Malthusian parameter or intrinsic 
growth rate. In other words, asynchronous exponential growth means that solutions tend 
to a fixed size distribution (often called the stable size profile) after a suitable rescaling 
of the semigroup. 

In this section we show that solutions of our model exhibit the same asymptotic 
behavior. This question can be addressed effectively in the framework of semigroup 
theory (see [U IH IS]). Briefly, to establish balanced exponential growth, one needs to 
show that the growth bound of the semigroup is governed by a leading eigenvalue of finite 
(algebraic) multiplicity of its generator, and there exists a spectral gap, i.e. the leading 
eigenvalue is isolated in the spectrum. Our first result will assure this latter condition. 
Moreover, if it is also possible to establish that the semigroup is irreducible then one has 
that the algebraic multiplicity of the spectral bound equals one (with a corresponding 
positive eigenvector) and the semigroup exhibits asynchronous exponential growth. 

Lemma 3.3. The spectrum of A can contain only isolated eigenvalues of finite algebraic 
multiplicity. 

Proof. We show that the resolvent operator R(X, A) is compact. Since A is a bounded 
perturbation of A it is enough to show that R(X,A) is compact. This follows however, 
from the regularity of the solution of the resolvent Equation (|2.10[) and noting that 
C 2 [0,m] C W^^m) © K 2 which is compactly embedded in L x (Q,m) © R 2 by the 
Rellich-Kondrachov Theorem [TJ Theorem 6.3, Part I]. The claim follows on the ground 
of [5J Proposition II.4.25]. □ 
Next we recall some necessary basic notions from linear semigroup theory, see e.g. 
EJ[6]. A strongly continuous semigroup {<S(t)} t>0 on a Banach space y with generator 
O and spectral bound 

s (O) = sup { Rc(A) : A e <r (O) } 
is said to exhibit balanced exponential growth if there exists a projection II on y such 
that 

lim \\e- s{0)t S{t) -n|| =0. 

t— >oo 

If the projection P is of rank one then the semigroup has asynchronous exponential 
growth. Moreover, the growth bound ujq is the infimum of all real numbers u> such that 
there exists a constant M > 1 with ||<S(t)|| < Me"'. We also recall (see e.g. [5J C- 
III Definition 3.1]) that a positive semigroup {<S(i)}t>o on a Banach lattice y is called 
irreducible if there is no S(t) invariant closed ideal of y except the trivial ones, {0} and 

y. 

Theorem 3.4. Model (ll.l[) - (|1.3l) admits a finite dimensional global attractor. 

Proof. We have shown in Theorem 12.11 that solutions of our model are governed by a 
positive semigroup. Derndingcr's Theorem (see e.g. [5J Theorem VI. 1.15]) implies that 
the spectral bound of the generator equals to the growth bound of the semigroup, i.e. 
s(A) = ujq. Lemma 3.1 implies that the spectral bound s(A) is an eigenvalue of finite 
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algebraic (hence geometric) multiplicity unless the spectrum is empty. If the spectrum of 
A is empty we have by definition ujq = — oo and every solution tends to zero. Otherwise 
we have for the semigroup {T(t)} t>0 generated by the closure of A 



lim 

t— f oo 



-s(A)t, 



T(t)-U = 0, (3.17) 



where II is the projection onto the finite dimensional eigenspace corresponding to the 
eigenvalue s(A). □ 
Recall that a subspace / C Z of a Banach lattice Z is an ideal iff / G / and \g\ < \ f\ 
implies that g £ I. The following useful result is from [2] C-III, Proposition 3.3]. 

Proposition 3.5. Let B be the generator of a positive semigroup lA(t) on the Banach 
lattice Z and K a bounded positive operator. Let V(t) be the semigroup generated by 
B + K . For a closed ideal L C Z the following assertions are equivalent: 

(i) L is V -invariant, 

(ii) I is invariant under both U and K . 

Wc introduce the recruitment operator K = A — A. 

Theorem 3.6. Assume that (|2.9[) holds true and (3 > 0. Then the semigroup generated 
by the closure of A exhibits asynchronous exponential growth. 

Proof. Our goal is to apply Proposition 13.51 for the operator A, whose closure is the 
generator of a positive semigroup as shown in Theorem 12.11 and to K, which is clearly 
positive and bounded. Every closed ideal I of X can be written as li ffi/2 ffi-^3, where I\ is 
a closed ideal in the Banach lattice L 1 (0, m) and I2, 13 are closed ideals in R. Note that, R 
admits only two ideals, i.e. {0} or R itself. Next we observe that non-trivial closed ideals 
in I/ 1 (0, to) can be characterized via closed subsets G of positive measure of (0, to). That 
is, the subspace J is a closed ideal of i 1 (0,m) if it contains the functions / 6 L 1 (0,to) 
vanishing on G. Next we show that no non-trivial closed ideal I = L\ © L2 is invariant 
under K or under the semigroup generated by A. If I\ 7^ {0} then the condition /3 > 
guarantees that Ku(s) = JJ™ 0(s,y)u(y)dy > 0, for every s <E (0,m) for any u € I\, i.e. 
the image Ku does not vanish anywhere, hence by the previous characterization we must 
have I\ = L 1 (0, to). Moreover, in this case j3 > implies that we must have I2 = I3 = R, 
since /J™ /3(0, y)u{y) dy > and J™ (3{m,y)u{y) dy > for any u ^ 0. On the other 
hand, if h = {0} then we have D(A) Hi = {0}, hence the restriction of A to R © R (or 
even to R) does not generate a semigroup. This means that / cannot be invariant under 
the semigroup generated by A. That is we have by Proposition 13.51 that the semigroup 
generated by A has no non-trivial closed invariant ideal. Therefore it is irreducible, and 
solutions exhibit asynchronous exponential growth, see e.g. [4]. □ 

4. Concluding remarks 



In this note we introduced a linear structured population model with diffusion in the 
size space. Introduction of a diffusion is natural in the biological context [19| . since unlike 
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in age-structured models, individuals that have the same size initially, may disperse as 
time progresses. In other words, diffusion amounts to adding noise in a deterministic 
fashion. We equipped our model with generalized Wentzcll-Robin boundary conditions. 
We showed that the model is governed by a positive quasicontractivc semigroup on the 
biologically relevant state space. Furthermore we have characterized the asymptotic 
behavior of solutions via balanced exponential growth of the governing semigroup. We 
also established that solutions exhibit asynchronous exponential growth if the function 
j3 is strictly positive. The question of irreducibility of the semigroup generated by the 
Laplace operator with mixed Robin boundary conditions on a L p -space (for 1 < p < oo) 
was addressed in the recent work by Haller-Dintclmann et al. [23]. It is expected that 
a similar result holds if generalized Wentzell-Robin boundary conditions are imposed. 
This is a topic of ongoing research. 

To our knowledge, generalized Wentzcll boundary conditions have not yet been intro- 
duced in the context of physiologically structured population dynamical models. Their 
power is to model sources and sinks on the boundary of a domain [16j which normally is a 
set of measure zero and therefore seems to play no role in an integral term. Interestingly, 
sinks on the boundary can cause ill-posedncss in space dimensions > 2 as Vazquez and 
Vitillaro [22 have shown. 

In the future we will extend our model to incorporate interaction variables, to allows 
competition. Then model (jl.ip ~ p.3p becomes a nonlinear one, and the mathematical 
analysis will become more difficult. To our knowledge, positivity results are rather rare 
in the literature for nonlinear models. In |llj it was shown that the nonlinear semigroup 
generator satisfies the positive minimum principle, hence the semigroup is positive. This 
however, does not apply to population models, as the positive cone of the natural state 
space L 1 has empty interior, hence the positive minimum principle docs not apply. 
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Abstract. We consider a linear size-structured population model with diffu- 
sion in the size-space. Individuals are recruited into the population at arbi- 
trary sizes. We equip the model with generalized Wentzcll-Robin (or dynamic) 
boundary conditions. This approach allows the modelling of populations in 
which individuals may have distinguished physiological states. We establish 
existence and positivity of solutions by showing that solutions are governed 
by a positive quasicontractive semigroup of linear operators on the biologically 
relevant state space. These results are obtained by establishing dissipativity of 
a suitably perturbed semigroup generator. We also show that solutions of the 
model exhibit balanced exponential growth, that is, our model admits a finite- 
dimensional global attractor. In case of strictly positive fertility we arc able to 
establish that solutions in fact exhibit asynchronous exponential growth. 

1. Introduction. A significant amount of interest has been devoted to the analysis 
of mathematical models arising in structured population dynamics (see e.g. [21, 27] 
for references). Such models often assume spatial homogeneity of the population 
in a given habitat and only focus on the dynamics of the population arising from 
differences between individuals with respect to some physiological structure. In this 
context, reproduction, death and growth characterize individual behavior which 
may be affected by competition, for example for available resources. 

In a recent paper, Hadeler [18] introduced size-structured population models 
with diffusion in the size-space. The biological motivation is that diffusion allows 
for "stochastic noise" to be incorporated in the model equations in a deterministic 
fashion. The main question addressed in [18] is what type of boundary conditions 
are necessary for a biologically plausible and mathematically sound model. In this 
context some special cases of a general Robin boundary condition were consid- 
ered. Diffusion terms have been introduced into structured population models by 
Waldstattcr et al. [26], Milncr and Patton [22] and Langlais and Milncr [20] in the 
context of host-parasite models, where the parasite load is the continuous structure 
variable. For example, when the structuring variable represents a parasite load as 

2000 Mathematics Subject Classification. 92D25, 47N60, 47D06, 35B35. 

Key words and phrases. Structured populations; diffusion; Wentzcll-Robin boundary condi- 
tion; semigroups of linear operators; spectral methods; stability. 
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in [26], at least one special compartment of individuals arises, namely the class of 
the uninfected ones. In a model where the structuring variable is continuous this 
state corresponds to a set of measure zero, hence it does not carry mass. The true 
meaning and advantage of employing Wcntzell boundary conditions is that it allows 
this special state to carry mass. 

In this paper we introduce the following linear size-structured population model 

u t (s,t) + (j(s)u(s,t)) 8 

= (d(s)u s (s,t)) s - fi(s)u(s,t)+ p(s,y)u(y,t)dy, se(0,m), (1) 

Jo 

[(d(s)u s (s,t)) s ] s=0 - b u s (0,t) +c u(0,t) = 0, (2) 

[(d(s)u s (s,t)) s ] s=m + b m u s (m,t) + c m u(m,t) = 0, (3) 

with a suitable initial condition. The function u = u(s,t) denotes the density of 
individuals of size, or other developmental stage, s at time t. Note that, we use 
as the minimal value of the variable of s only for mathematical convenience. It 
may well be replaced by an arbitrary value s m i n . The non-local integral term in 
(1) represents the recruitment of individuals into the population. Individuals may 
have different sizes at birth and /3(s, y) denotes the rate at which individuals of 
size y "produce" individuals of size s. Further biologically relevant assumptions 
may be made on the fertility, such as fi{s,y) = for s > y, i.e. parents cannot 
have larger offspring. We also note that from this general model a single state at 
birth model may be deduced by formally replacing the fertility function j3 with an 
appropriate delta function, see e.g. [21]; Chapter I Section 4. fi denotes the size- 
specific mortality rate while 7 denotes the growth rate, d stands for the size-specific 
diffusion coefficient, which we assume to be strictly positive. 60 an d b m are positive 
numbers, while cq and c m are non-negative. We will discuss special values for these 
constants in Equation (5) below. 

Equation (1) describes the evolution of a "proper" size-structured population, in 
contrast to the one where it is assumed that all newborns enter the population at 
a minimal size. In that case, assuming that the growth rate is positive, i.e. individ- 
uals do not shrink, the linear model can be rewritten as an age-structured model, 
see e.g. [21]. We also refer the interested reader to [3, 7, 10] where size-structured 
models with distributed recruitment processes were investigated. We also note 
that although we refer to the structuring variable s as size, it could well represent 
any other physiological characteristic of individuals such as accumulated energy or 
biomass, volume etc. The boundary conditions (2)-(3) are the so called general- 
ized Wentzell-Robin or dynamic boundary conditions. These "unusual" boundary 
conditions were investigated recently for models describing physical processes such 
as diffusion and wave propagation, see e.g. [11, 12, 16]. Briefly, they are used to 
model processes where particles reaching the boundary of a domain can be either 
reflected from the boundary or they can be absorbed. Hence the boundary points 
can carry mass. Our goal here is to introduce this rather general type of boundary 
condition in the context of models describing the evolution of biological populations, 
with particular focus on positivity and asymptotic behavior of solutions. Potential 
applications include cell populations with resting states at s = and s = m, or 
models for populations structured by an infection level or parasite load [26]. 

The first works introducing boundary conditions that involve second order deriva- 
tives for parabolic or elliptic differential operators go back to the 1950s, see the 
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papers by Feller [13, 14] and Wcntzell [24, 25]. These first studies were purely mo- 
tivated from the mathematical point of view. The original question was to identify 
the set of all possible boundary conditions that give rise for a parabolic differential 
operator to generate a contraction semigroup on an appropriate state space. The 
abstract mathematical analysis gave a clue where physical intuition failed previ- 
ously. 

The boundary conditions (2)-(3) are in general form, and we shall now specify 
the constants bo,b m and co,c m to give a biological explanation for the boundary 
conditions. Integration of Equation (1) from to m yields for U(t) = J Q u(s,t)ds 

^ £7(t) = 7(0)it(0,t) - j(m)u(m,t) + d(m)u s (m,t) - d(0)u s (0,t) 

+ / y)u(y,t)dyds - / fx(s)u(s,t)ds 

Jo Jo Jo 

=: B(t) - D(t) + 7(0)u(0, t) - j(m)u(m,t) + d(m)u s (m,t) - d(0)u s (0,t), 

where B and D denote the combined birth and death processes, respectively We 
note that formally, by replacing the diffusion term by its counterpart from Equation 
(f ), the boundary conditions (2)- (3) can be cast in the dynamic form 



u t (0, t) =u(0, t)(-7'(0) - MO) - c ) + u s (0, t)(bo - 7(0)) + / P{0, y)u(y, t) dy, 

Jo 

u t (m,t) =u(m,t)(-^' (m) - /i(m) - c m ) + u s (m,t) (—b m - 7(m)) 

+ / P{m 1 y)u{y,t)dy. 
Jo 

(4) 

These are the governing equations for individuals of minimum and maximum sizes, 
respectively. It is natural to assume that in the absence of mortality and recruit- 
ment, i.e. when B(-) = D(-) = 0, the total population size U(t) + u(0,t) + u(m,t) 
remains constant at every time t. Mathematically, this amounts to the condition 

0=— U(t)+u t (0,t)+u t (m,t) 

=u(0, t) ( 7 (0) - 7 '(0) - c ) + u{m, t) (-7(m) - c m - 7 '(m)) 

+ u s (0, t) (b - d(0) - 7 (0)) + u s (m, t) (d{m) - 7 (m) - b m ) . 

In order to guarantee conservation of total population in the absence of birth and 
death processes, we make the following assumptions 

co = 7(0) - Y(0), Cm = - 7 (m) - Y(m), b = d(0) + 7 (0), b m = d(m) - 7 (m). (5) 

We note that condition (5) together with the assumption that b is positive and c is 
non-negative impose a restriction on the growth rate 7. 

The dynamic boundary conditions (4) can now be written as 



Ut(0,*) = (-A*(0)-7(0))«(0,*)+d(0K(0,t)+ / P(0,yMy,t)dy, 

J ° fm (6) 
ut(m,t) = (7(m) — ii(m))u(m,t) — d(m)u s (m,t) + / (3{m,y)u(y,t)dy. 



Hence the dynamics of individuals in the two special states s = and s = m are 
governed by equations (6). The meaning of the governing equations (6) is intuitively 
clear. For example, at s = 0, individuals are leaving this compartment due to growth 
and mortality and are recruited according to the integral term, while the diffusion 
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accounts for the flux through this state (this also yields a loss term, since the outer 
normal derivative is — u s ). 

We impose the following assumptions on the model ingredients 

fiEC([0,m]), p £ C([0,m] x [0,m]), /3, p, > 0, 7, d £ C x ([0, m]), d > 0. 

In this note first we establish existence and positivity of solutions of model (1)- 
(3). This existence proof follows similar arguments developed in [11, 12]. The 
significant difference is that solutions to our model are not necessarily governed 
by a contraction semigroup, hence as in [8, 9] we need to rescale the semigroup to 
obtain the dissipativity estimate. As a result of the dissipativity calculation, we 
show that the resolvent operator of the semigroup generator is positive. This was 
not established in [12]. 

In Section 3 we investigate the asymptotic behavior of solutions. First we es- 
tablish that solutions of the model equations exhibit balanced exponential growth, 
in general. This is an interesting phenomenon, often observed for linear structured 
population models, see e.g. [8, 10]. In some sense, it is a stability result and char- 
acterizes the global asymptotic behavior of solutions to the model. Then, assuming 
that fertility is strictly positive, we are able to show that after a rescaling by and 
exponential factor solutions actually tend to a fixed size-distribution. This is shown 
via establishing irrcducibility of the governing semigroup. 

2. Existence and positivity of solutions. In this section we are going to es- 
tablish the existence of a positive quasicontractive semigroup of operators which 
governs the evolution of solutions of (l)-(3). For basic definitions and results used 
throughout this paper we refer the reader to [1, 2, 4, 6]. Let 

X = (1^(0, m) ®M 2 , || 
where for (x, xq, x m ) £ X the norm is given by 

||(x,x ,a; m )||^ = \\x\\ L i +ci\x \ + c 2 |a; m |, (7) 

for some positive constants c\ and c 2 that we will specify later. Then X is a 
Banach lattice. We identify a function u £ C[0, m] with its restriction triple 
(«|(o, m ),u(0),«(m)) £ X. With this identification, the set C 2 [0,to] is dense in 
X with respect to the || ■ H^-norm. Let 

D[A) = {u £ C 2 [0,m] : *u £ L\0, m), {d(s)u'(s))'\ a=Q - b o u'(0) + cqu(0) = 0, 
(d(s)u'(s))'\ s=m + b m u'(m) + c m u(m) = 0} , 

where 

Vu(s) = (d(s)u'(s))' - (7(s)u(s))' - fJ.(s)u(s) + / /3(s, y)u(y) dy. 



The operator A with domain D(A) is then defined by 
/ ^>u 

Au = (b - 7 (0))£u | s=0 + / m /3(0, y)u(y) dy - Po u(0) 

\(-b m - 7N)|« \ s=m + J™ f3(m, y)u(y) dy - p m u(m)j 
where we have set 

Po = /x(0) + c + V(0), p m = fi(m) + Cm + j'{m), 
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for short. We use [6, Theorem 3.15] (see also [2, Section A-II, Theorem 2.11]) which 
characterizes generators of contractive semigroups via dissipativity. Wc establish 
the existence of a quasicontractivc semigroup for the general boundary condition (2)- 
(3). Recall that a strongly continuous semigroup {T(t)} t>0 is called quasicontr active 
if ||T(i)|| < e ut for some uj G K and it is called contractive if the choice cj < is 
possible. A linear operator A with domain D(A) is dissipative, if for all x G D(A) 
and A > one has |(AZ — A)x\ \ > A||x||. 

Theorem 2.1. Assume that 

d(0) , dim) 

a = 7 1 > 0, and c 2 = \> >0. 8 
oq - 7(0) J(m) + b m 

Then the closure of the operator A is the infinitesimal generator of a positive qua- 
sicontractive semigroup of bounded linear operators on the state space X ( where the 
weights in equation (7) are chosen accordingly). 

Proof. We introduce the modified operator A with j3 — in the definition of A. For 
A > and h ^ X,u <E D(A) we consider the equation 

u - A (A - ljT) u = h, on [0,m]. (9) 

That is 

h(s) = u(s) — A ((d(s)u'(s))' — (j(s)u(s))' — /j,(s)u(s) — uju(s)j , s G (0, m), 

(10) 

A" 1 /^) = u'(0)(7(0) - &o) + "(0) (A" 1 + 7'(0) + m(0) + c + w) , (11) 
\~ l h{ra) = u'(m)(7(m) + b m ) + u(m) (A -1 + j'(m) + fi(m) + c m + uj) . (12) 

Next we multiply Equation (10) by Xu+( s )i where \u+ denotes the characteristic 
function of u + , and integrate from to m. The boundary d[u > 0] consists of two 
parts, Ti = d\u > 0] (~l (0,m) and r2 = [u > 0] fl {0,m}. The term d(s)u'(s) gives 
negative contributions on Ti since the outer normal derivative of it is negative. The 
term j(s)u(s) gives no contributions on Ti since u = there. Hence we obtain 

||w + ||i <A (8gn(u + (m))d{m)u'(m) - sgn(u+(0))d(0)u'(0)) 
— A (sgn(ii + (m))7(m)w.(m) — sgn(u + (0))7(0)u(0)) 



-A/ Xu+(s)K s ) u ( s )& s - Aw||it + ||i + / \u+(s)H s )d s - 
Jo Jo 

This, combined with equations (11)-(12), yields 

||u + ||i + sgn(ii + (m))u(m) [ ^ m — (1 + A h'(m) + /u(m) + c m + w)) + A7(m) ] 

\"f(m) + b m ) 

+ sgn( W +(0)) U (0) ( 6o d( ° } (0) (1 + A ( 7 '(0) + M (0) + co + w)) - A 7 (0)^ 
< -A / Xu+(s)/u(s)u(s)ds - Aw||u + ||i + / y„+ (s)fr(s) ds 



+ sgn(u+(m))/i(m) + sgn(u+(0))/i(0)- 



7(m) + 6 m 6 - 7(0) 

(13) 
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Similarly, multiplying Equation (10) by — \ u - (s) and integrating from to to, we 
obtain 

- sgn(u~ (m))u(m) ( (1 + A (Y(m) + /u(m) + c m + u)) + Xj(m) ) 

\7{m) + b m ) 

- sgn( M - (0))u(0) ( bQ d }°\ 0) (1 + A (7'(0) + MO) + c + w)) - A 7 (0)^ 

/> m /> rn 

<A/ x«-( s )M s M s ) ds - Aw||u"||i - / Xu-(s)h(s)ds 



sgn(u-(m))h(m) , T\ - sgn(«-(0))ft(0)- d(0) 



(14) 



7(to) + 6 m & - 7(0) ' 

Adding (13) and (14) yields 

Nil + l«(0)| L o d }^ {0) C 1 + A (7'(0) + M(0) + co + w)) - A 7 (0)) 
/ dim) 

+ l u ( TO )l ? I ri (! + A (V( m ) + MM + c '» + w )) + A 7(™) 
\7(m) + o m 

(X U+ 00 ~ Xu- (*)) /»(«) ds + ft(0) , M« + (0)) - sgn(u-(0))) (15) 

oo - 7l u J 

d( JTL ) 

- Aw||u||i + fe(m) . (sgn(u+(m)) - sgn(u~(m))) 

/•7TJ 

-X (x«+(s) - X«-(s))M( s )"(s)ds. 
Jo 

Assuming that condition (8) holds true the left hand side of inequality (15) can be 
estimated below, for w £ K large enough, by 



< 



/r ,M d (0) 1 I / \ I g?(to) 

\u\\i + |u(0)|- —r + \u[m)\- 



b Q - 7(0) 7(771) + b. m ' 

Similarly, the right hand side of inequality (15) can be estimated above by 

INIi + IMO)lr^W + K-)l- d{m) 



6o-7(0) l{m) + b m 

Hence, if condition (8) is satisfied we have the dissipativity estimate 



Ml* < \\h\\x = 



-XiA-wI 



X 



for the operator A — loI. 

For the range condition we need to show that whenever h is in a dense subset of 
X then the solution u of Equation (9) belongs to the domain of A. Since we assumed 
that d > 0, i.e. we have true Wentzell boundary conditions at both endpoints of the 
domain, the required regularity of the solution u £ C 2 [0,to] and hence the range 
condition follows from [15, Theorem 6.31]. Thus the closure of A — ulL is a generator 
of a contractive semigroup by [6, Theorem 3.15]. 

Next we observe that for h £ X + every term on the right hand side of Equation 
(14) is non-positive, while every term on the left hand side of Equation (14) is non- 
negative. The inequality can only hold for h £ X + if = u(0) = u(m) = 0. 
This proves that the resolvent operator R(X, A — ljI) = (XT— (A — uil))^ 1 is positive 
for A large enough, hence the closure of A — uT generates a positive semigroup. A 
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simple perturbation result yields that the closure of A is a generator of a positive 
quasicontractive semigroup. 

Finally we note that, the operator A — A is positive, bounded and linear, hence 
it generates a positive semigroup {S(t)} t>0 which satisfies 



The proof of the Theorem is now completed on the grounds of the Trotter product 



Remark 1. Conditions (8) are clearly satisfied if only the denominators are posi- 
tive. However, the notation in (8) gives immediately the norm on the state space. 

Remark 2. We note that the biologically natural assumptions in (5) that guaran- 
tee conservation of the total population in the absence of mortality and recruitment 
imply that c\ = ox = 1. Hence the mathematical calculations coincide with the bio- 
logical intuition. It can be shown that in the absence of mortality and recruitment 
the operator A generates a contraction semigroup on the state space (X, || • \\x), in 
fact ||T(i)|U = 1 for all t > 0. 

3. Asymptotic behavior. In this section we investigate in the framework of semi- 
group theory the asymptotic behavior of solutions of model (l)-(3). Since our model 
is a linear one, we may expect that solutions either grow or decay exponentially, 
unless they arc at a steady state. For simple structured population models in fact it 
is often observed (see e.g. [5, 8, 10, 17]), that solutions grow exponentially and tend 
to a finite-dimensional (resp. one-dimensional) global attractor. This phenomenon 
is called balanced (resp. asynchronous) exponential growth. The rate of exponen- 
tial growth is called the Malthusian parameter or intrinsic growth rate. In other 
words, asynchronous exponential growth means that solutions tend to a fixed size 
distribution (often called the stable size profile) after a suitable rescaling of the 
semigroup. 

In this section we show that solutions of our model exhibit the same asymptotic 
behavior. This question can be addressed effectively in the framework of semigroup 
theory (see [2, 4, 6]). Briefly, to establish balanced exponential growth, one needs to 
show that the growth bound of the semigroup is governed by a leading eigenvalue 
of finite (algebraic) multiplicity of its generator, and there exists a spectral gap, 
i.e. the leading eigenvalue is isolated in the spectrum. Our first result will assure 
this latter condition. Moreover, if it is also possible to establish that the semigroup 
is irreducible then one has that the algebraic multiplicity of the spectral bound 
equals one (with a corresponding positive eigenvector) and the semigroup exhibits 
asynchronous exponential growth. 

Lemma 3.1. The spectrum of A can contain only isolated eigenvalues of finite 
algebraic multiplicity. 

Proof. We show that the resolvent operator i?(A, ^4) is compact. Since A is a 
bounded perturbation of A it is enough to show that R(X,A) is compact. This 
follows however, from the regularity of the solution of the resolvent Equation (9) and 
noting that C 2 [0, m] C W 1 ' 1 (0, m) ffiR 2 which is compactly embedded in L 1 (0, m) © 
M 2 by the Rcllich-Kondrachov Theorem [1, Theorem 6.3, Part I]. The claim follows 
on the ground of [6, Proposition II. 4. 25]. □ 



l|5(*)||<e : 



IB 



where B = 



OO * 




□ 
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Next we recall some necessary basic notions from linear semigroup theory, see 
e.g. [2, 4, 6]. A strongly continuous semigroup {S(t)} t>0 on a Banach space y with 
generator O and spectral bound 

s (O) = sup { Rc(A) : A G a (O) } 

is said to exhibit balanced exponential growth if there exists a projection II on y 
such that 

lim \\e- s{0)t S{t) -n|| =0. 

t— >oo 

If the projection P is of rank one then the semigroup has asynchronous exponential 
growth. Moreover, the growth bound loq is the infimum of all real numbers u> such 
that there exists a constant M > 1 with ||<S(t)|| < Me" 1 . We also recall (see e.g. [2, 
C-III Definition 3.1]) that a positive semigroup {S(t)} t>0 on a Banach lattice y is 
called irreducible if there is no S(t) invariant closed ideal of y except the trivial 
ones, {0} and y. 

Theorem 3.2. Model (l)-(3) admits a finite dimensional global attractor. 

Proof. We have shown in Theorem 2.1 that solutions of our model are governed 
by a positive semigroup. Dcrndingcr's Theorem (see e.g. [6, Theorem VI. 1.15]) 
implies that the spectral bound of the generator equals to the growth bound of the 
semigroup, i.e. s(A) = ujq. Lemma 3.1 implies that the spectral bound s(A) is an 
eigenvalue of finite algebraic (hence geometric) multiplicity unless the spectrum is 
empty. If the spectrum of A is empty we have by definition uq = — oo and every 
solution tends to zero. Otherwise we have for the semigroup {T(t)} t>0 generated 
by the closure of A 

lim e' s{A)t r{t) n =0, (16) 

t— > oo 

where LT is the projection onto the finite dimensional eigenspace corresponding to 
the eigenvalue s(A). □ 
Recall that a subspace / C Z of a Banach lattice Z is an ideal iff / G / and 
\g\ < |/| implies that g G I. The following useful result is from [2, C-III, Proposition 
3.3]. 

Proposition 3.3. Let B be the generator of a positive semigroup U(t) on the Ba- 
nach lattice Z and K a bounded positive operator. Let V(t) be the semigroup gen- 
erated by B + K . For a closed ideal L C Z the following assertions are equivalent: 

(i) I is V -invariant, 

(ii) L is invariant under both hi and K . 

We introduce the recruitment operator K = A — A. 

Theorem 3.4. Assume that (8) holds true and {3 > 0. Then the semigroup gener- 
ated by the closure of A exhibits asynchronous exponential growth. 

Proof. Our goal is to apply Proposition 3.3 for the operator A, whose closure is the 
generator of a positive semigroup as shown in Theorem 2.1 and to K , which is clearly 
positive and bounded. Every closed ideal L oi X can be written as I\ ®L 2 ®L' il where 
L\ is a closed ideal in the Banach lattice L 1 (0,to) and l2,Iz are closed ideals in R. 
Note that, R admits only two ideals, i.e. {0} or R itself. Next we observe that non- 
trivial closed ideals in L 1 (0, m) can be characterized via closed subsets G of positive 
measure of (0, to). That is, the subspace J is a closed ideal of L 1 (0, m) if it contains 
the functions / G -^(O, to) vanishing on G. Next we show that no non-trivial closed 
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ideal I = I\ (Bh ®h is invariant under K or under the semigroup generated by A. If 
I\ ^ {0} then the condition (3 > guarantees that Ku(s) ~ f Q (3(s, y)u(y) dy > 0, 
for every s £ (0,m) for any u £ 7i, i.e. the image Ku does not vanish anywhere, 
hence by the previous characterization we must have I\ = ^(O^m). Moreover, in 
this case (3 > implies that we must have I2 = I3 = R, since J j3(0, y)u(y) dy > 
and /J™ /3(m, y)u(y) dy > for any u ^ 0. On the other hand, if I\ = {0} then we 

have D(A) D I = {0}, hence the restriction of A to R © R (or even to R) does not 
generate a semigroup. This means that I cannot be invariant under the semigroup 
generated by A. That is we have by Proposition 3.3 that the semigroup generated by 
A has no non-trivial closed invariant ideal. Therefore it is irreducible, and solutions 
exhibit asynchronous exponential growth, see e.g. [4]. □ 
The previous result completely characterizes the asymptotic behavior of solutions 
to the population model. That is, solutions behave asymptotically as e rt u*(s) in- 
dependently of the initial condition, where r is the so called Malthusian parameter, 
and is often referred to as the final size distribution. 



4. Concluding remarks. In this note we introduced a linear structured popula- 
tion model with diffusion in the size space. Introduction of a diffusion is natural in 
the biological context [18], since unlike in age-structured models, individuals that 
have the same size initially, may disperse as time progresses. In other words, dif- 
fusion amounts to adding noise in a deterministic fashion. We equipped our model 
with generalized Wcntzell-Robin boundary conditions. We showed that the model is 
governed by a positive quasicontractive semigroup on the biologically relevant state 
space. Furthermore we have characterized the asymptotic behavior of solutions 
via balanced exponential growth of the governing semigroup. We also established 
that solutions exhibit asynchronous exponential growth if the function /3 is strictly 
positive. An important biological consequence of asynchronous exponential growth 
is population stabilization in the sense that the proportion of the population in 
any subset of the structure space converges to a limiting value as time evolves, 
independently of the initial state of the population. The question of irreducibility 
of the semigroup generated by the Laplace operator with mixed Robin boundary 
conditions on a L p -spacc (for 1 < p < 00) was addressed in the recent work by 
Hallcr-Dintclmann et al. [19] . It is expected that a similar result holds if general- 
ized Wcntzell-Robin boundary conditions arc imposed. This is a topic of ongoing 
research. 

In our model we have taken the view that individuals may be recruited into the 
population at different sizes. This appears to be the natural choice in the context of 
general physiologically structured population models, as opposed to age-structured 
models, where every individual is born at the same age zero. It is interesting to 
investigate whether a mathematically sound "limiting relationship" exists between 
models with infinite states at birth and one state at birth. This will be addressed 
in future work. 

The power of generalized Wentzcll boundary conditions in the context of pop- 
ulation models is to allow the boundary states to carry mass. This is especially 
interesting in the L p context as the boundary is a set of measure zero and therefore 
seems to play no role in an integral term. Interestingly, sinks on the boundary 
can cause ill-posedness in space dimensions > 2 as Vazquez and Vitillaro [23] have 
shown. 
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In the future we will extend our model to incorporate interaction variables, to 
allow competition. Then model (l)-(3) becomes a nonlinear one, and the math- 
ematical analysis will become more difficult. To our knowledge, positivity results 
are rather rare in the literature for nonlinear models. In [11] it was shown that the 
nonlinear semigroup generator satisfies the positive minimum principle, hence the 
semigroup is positive. This however, does not apply to population models, as the 
positive cone of the natural state space L 1 has empty interior, hence the positive 
minimum principle does not apply. It will be also interesting to consider more gen- 
eral models with a finite number of structuring variables, such as age-size structured 
models. Then the domain will be a cube [0, 1]" and the prescription of appropriate 
boundary conditions will be much more involved. 
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